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Recently it has been discovered that lanthanum, yttrium, and other metal hydride films show 
dramatic changes in the optical properties at the metal-insulator transition. Such changes on a high 
energy scale suggest the electronic structure is best described by a local model based on negatively 
charged hydrogen (H~) ions. We develop a many-body theory for the strong correlation in a H - 
ion lattice. The metal hydride is described by a large (7-limit of an Anderson lattice model. We 
use lanthanum hydride as a prototype of these compounds, and find LaH3 is an insulator with 
a substantial gap consistent with experiments. It may be viewed either as a Kondo insulator or 
a band insulator due to strong electron correlation. A H vacancy state in LaH3 is found to be 
highly localized due to the strong bonding between the electron orbitals of hydrogen and metal 
atoms. Unlike the impurity states in the usual semiconductors, there is only weak internal optical 
transitions within the vacancy. The metal-insulator transition takes place in a band of these vacancy 
states. 

PACS numbers: 71.30.+h, 71.15.Mb, 71.55.Ht, 72.15.Eb 



I. INTRODUCTION 

Huiberts et al. |lj have recently reported dramatic 
changes in the optical properties of lanthanum, yttrium, 
and other rare earth hydride films Q with changing hy- 
drogen content. This phenomenon has also been recently 
observed in other metal hydrides, such as Gd-Mg alloy 
hydrides ||. By changing the hydrogen gas pressure, or 
by electrochemical means || the films can be continu- 
ously and reversibly switched from a shiny mirror (good 
metal) to a transparent window (insulator) in a fraction 
of a second. Although many metal-insulator transitions 
are known, this type of switchable optical phenomena is 
very unusual, and potentially of considerable technolog- 
ical importance since the transition leads to spectacu- 
lar effects in the visible light region. For example, this 
class of materials is very different from the previously re- 
ported charge density wave compound KM0O3. In that 
compound, the frequency dependent conductivity may be 
changed by replacing some of Mo atoms by W but 
the change in reflectivity is in the infrared or lower fre- 
quency region, while the metal hydrides show dramatic 
changes in the visible. Since the hydrides are easily tun- 
able, this class of compounds is also ideal for basic re- 
search to understand the metal-insulator transition in 
general. Furthermore, since the optical switching is re- 
alized at room temperature and at normal pressure, and 
since there appears to be considerable scope for short- 
ening the time-scale of the transition through chemical 
or electrochemical means Q , the phenomenon is very at- 
tractive for optical devices. 

As Huiberts et al. jjj reported, the transitions mea- 
sured in the optical transmission and in the electric 



resistivity appear at the same hydrogen concentration 
x = x c ~ 2.80 for both the YH X and LaH x . The optical 
gap in the insulating phase is ~ 1.8 eV for LaH3 and 
~ 2.8 eV for YH3. The transition is clearly of electronic 
origin in the lanthanum hydrides where the crystal struc- 
ture remains face centered cubic (fee) for 2 < x < 3. In 
the yttrium hydrides, there is evidence that the observed 
metal-insulator transition is also of electronic origin || . 

The electronic structure underlying this behavior was 
poorly understood — indeed the standard local density 
approximation (LDA) calculations failed to predict a 
metal-insulator transition at all. In 1970, Switendick |l2j | 
used a non-self-consistent approximation to calculate the 
electronic structure of the YH3 , and found an energy gap 
of 1.5 eV. But more sophisticated state-of-the-art self- 
consistent LDA calculations by Dekker et al. (r|, Wang 
and Chou [[li} predict that LaH3 and YH3 are metals or 
semimetals. 

The structure changes are especially small in the LaH^ 
films, which may be considered as a prototype for this 
class of compounds. The La atoms always form a fee 
lattice with two H atoms occupying tetrahedrally coor- 
dinated sites. As x changes from 2 to 3, the octahedrally 
coordinated sites go from empty to fully occupied and 
a good metal evolves to a transparent insulator. The 
crystal structure is shown in Fig. |l|. The actual lattice 
constants of the dihydride and the trihydride are 10.7051 
a.u. (1 a.u.= 0.528 A) and 10.5946 a.u. respectively. The 
slight contraction due to the addition of H atom in LaH3 
indicates strong hybridization between La and H atoms. 
But the small difference in lattice constant is insignifi- 
cant, and will be neglected in our further consideration. 

In the metal hydrides, the hydrogen atom attracts one 
more electron and forms a negative hydrogen ion H _ . 



1 



The rare earth metal such as La has a formal valence 
3+. In a dihydride such as LaH2, the ls-orbital of the hy- 
drogen bands (H bands) are filled, leaving 1 electron per 
unit cell in the La conduction bands. This gives metallic 
behavior. In a trihydride such as LaH3, the H bands can 
hold all six valence electrons per unit cell. However, in 
the LDA calculations there is an overlap between the H 
and La 5e?-orbital bands, leading to a metal or semimetal 
rather than the observed transparent insulator. The H~ 
ion is a difficult case for the LDA, and a careful treatment 
of the correlation between the two electrons is required in 
order to obtain the bound H ion state with binding en- 
ergy of 0.7 eV (see next section). The poor treatment of 
correlations is most likely the reason why LDA is unable 
to explain the transparent insulating behavior of LaH3. 
It is well known that the LDA underestimates the en- 
ergy gap in semiconductors. This problem seems much 
more pronounced in rare earth hydrides: LDA predicts a 
metallic state for an insulator with substantial gap. 

There has been considerable theoretical activity since 
the original discovery ]l6| , p7| , pT|Jl9H2Ct | . In an early Letter 
p6| , we examined the role of electron correlation and pro- 
posed a local model to describe the electronic structure 
and the metal-insulator transition in these hydrides. The 
purpose of the present paper is to refine our calculations 
and to present a more complete theory for these hydrides. 
Recently, Eder et al. Jl?]] have proposed a theory based on 
the occupation-dependent hoppings on H sites and exam- 
ined a Kondo- like-insulator model for YH3 . Their theory 
has much in common with ours. Very recently Chang et 
al. [^0| reported GW method calculations for LaH3 and 
found a gap of ~ 0.5 eV. This again indicates the impor- 
tance of electron correlation. On the other hand, Kelly et 
al. JlSfl proposed that the insulating nature of YH3 may 
be explained within the LDA due to a more complicated 
hexagonal structure for Y. Their proposal, however, has 
difficulty explaining the insulating nature of many other 
hydrides, such as La hydrides, whose crystal structure is 
a always simple fee. 

In this paper, we examine the importance of electron 
correlation in metal hydrides, and develop a many-body 
theoretic framework to study their electronic structure. 
LaH3 may be viewed as a Kondo insulator if we 

start with metallic phase LaH2 and consider every addi- 
tional H atom to a be Kondo impurity. LaH3 may also 
be viewed as a band insulator due to the strong electron 
correlation, which suppresses the overlaps between the H 
and metal bands. When the H concentration is reduced 
in the LaH3 structure, the elementary entities are the 
H-vacancies. We find highly localized electronic states 
centered on these vacancies. The metal-insulator transi- 
tion takes place in a band of these vacancy states. The 
paper is organized as follows. 

In section II, we first review an appropriate wavefunc- 
tion for a single H~ ion in free space incorporating the 
electron correlation explicitly. We then extend the dis- 
cussion to a lattice of H ions, and use a microscopic 
model to calculate the electron motion (or the effective 



hopping integrals) in such a lattice. 

In section III, we propose that the metal hydrides are 
best described by a large [/-limit Anderson lattice Hamil- 
tonian, with U the electron Coulomb repulsion on the H 
atom. At every H site, electron state is either singly 
or doubly occupied. The electron hopping integrals be- 
tween H sites are described by those obtained in our 
microscopic calculations in section II. Other parameters 
are extracted from the LDA calculations. We use the 
Gutzwiller method to study this model Hamiltonian for 
the hydrides within this method the local constraint on 
the H sites is replaced by a set of renormalization factors. 

In section IV we present our results for the electronic 
structure of LaH2 and LaH3. The former is found to be 
a metal, and the latter an insulator, in agreement with 
experiments. 

Section V consists of the study of the localized H va- 
cancy state in LaH3 and the H impurity state in LaH2. 
The localized vacancy state is proposed to be the elemen- 
tary entity in the hydrides and to be responsible for the 
metal-insulator transition. We will use finite size calcula- 
tions to verify the symmetry and localized nature of the 
vacancy states. In section VI we use the calculated elec- 
tronic band structure to study the optical conductivity 
and density of states. We perform finite size calculations 
and compare the results to experimental data. 

II. ELECTRON CORRELATION IN 
NEGATIVELY CHARGE HYDROGEN IONS 

In this section we start by reviewing the electron cor- 
relation in a single negatively charged hydrogen ion (H _ 
ion) and the importance of the cooperative motion of the 
two electrons around the proton. We will then describe 
a local model for studying the many-body problem of a 
lattice of H~ ions. A microscopic model calculation will 
be presented to estimate the effective kinetic energy for 
the electrons in a H ion lattice. 



A. A single H ion in free space 

This is a venerable but interesting problem [pT]-p5|. 
It was carefully studied by astrophysicists back in the 
1950's, when the H - ion was found to be of great impor- 
tance for the opacity of the atmosphere of the sun and 
similar stars. Since hydrogen contains only one proton, 
the smallest charge of all nuclei, the Coulomb attrac- 
tion between an electron and its nucleus is the weakest 
and when an additional electron is added to a neutral H 
atom, the Coulomb repulsion between the two electrons 
becomes crucially important. Electron-electron correla- 
tion has the largest effect in H atoms compared to other 
atoms. As found by Chandrasakhar p2|| , Bethe and 
Salpeter | pl| , and Hylleraas et al. [231 ], a careful treat- 
ment of the correlation between the two electrons is nec- 
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essary to obtain the correct binding energy of about 0.7 
eV. This shows that the II ion is a bound state due to 
a strong electron correlation effect. 

Variational trial wavefunctions have been applied to 
the H - ion for there is no exact analytic solution for the 
ion. Wavefunctions with up to 24 variation parameters 
p3| were proposed to give the best estimate of the bind- 
ing energy. On the other hand, there is a simple but ex- 
cellent trial wavefunction introduced by Chandrasakhar 
in the 1940's [g2[ to describe the electron correlation in 
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where (ri,r2) are the electron coordinates with respect 
to the proton, x is t ne spin singlet spinor, and the con- 
stants a — 1.075, b = 0.478, c = 0.312 in atomic units. 
As illustrated in Fig. ||, the wavefunction can roughly 
be visualized as an inner electron of radius ~ 1, and an 
outer electron of radius ~ 2 orbiting around the pro- 
ton. The electron correlation is described by the term 
c|ri — r2 1 , which tends to keep the two electrons apart. 
This is similar to the Laughlin's wavefunctions for the 
fractional quantum Hall states, where electrons in high 
magnetic fields tend to stay apart and form a quantum 
liquid 26 1. One should notice that the second electron 
does not reside in the same orbital as the first electron 
although they have opposite spins, nor does it occupy the 
2s orbital which has a very different radial wavefunction. 
In this wavefunction the electrons are in a spin singlet 
state while the spin triplet state is unbound. The wave- 
function of Eq. (|l|) gives a ground state energy of H ion 
very close to the best estimate of Hylleraas |2j| using 24 
variational parameters. Note that the choice of the form 
of the variational wavefunction is very subtle and other 
simple choices mostly give an unbound electron. The 
standard LDA calculations for a single II do not give a 
bound state p7[ . 

The simple wavefunction of Eq. ([!]) will be used to 
develop a many body theory for the H~ ion lattice. The 
more refined 24 parameters wavefunction may produce 
a better result, but the calculation of hopping integrals 
involving such a complicated wavefunction will be very 
difficult and time consuming. 



B. H ions in crystal 

Since a single H - ion is a bound state, the ionic picture 
should be valid in certain hydride crystals. A local de- 
scription is suggested by the recently observed insulating 
gap in LaH3 and YH3. This led us to examine the effect 
of correlations on the H band width. In this section 
we shall focus on H~ ions. The effect of rare earth ions 
through their hybridization with H _ ions is also signifi- 
cant, and will be discussed in a later section. The form 
of the H~ bands in a many body theory is determined by 



the spectrum of the states obtained by removing an elec- 
tron from a lattice of H~ ions. We will use a local model, 
namely an orthogonal tight-binding model to describe the 
hole motion. Then the energy spectrum is determined by 
the effective hopping integrals t(d) of an electron from a 
n~ ion to a neutral H atom at distance d. Below we 
shall estimate t(d) from a microscopic model, in which 
the correlated Chandrasakhar wavefunction Eq. (|l|) will 
be employed. 

We consider a single pair consisting of one H ion and 
one H atom separated by a distance d as illustrated in 
Fig. § This results in a H 2 ion, an ion with three elec- 
trons moving around two protons. Previous studies have 
shown that ion is a bound negative ion [Q for d > 3 
a.u., the range of interest here (the nearest H-H distance 
separation in LaH3 is 4.58 a.u.). For d < 3 a.u., 
ion becomes unstable and dissociates into H2 and an un- 
bound, free electron. When a H~ ion and H atom are 
brought together from infinity, the ground state mani- 
fold of Hj" splits into odd(— ) and even(+) parity states 
with respect to the center of mass of the two protons. 
It is this splitting that determines the effective hopping 
matrix element t(d). Let E- and E+ be the energies of 
the odd and even parity states respectively, then we have 



t(d) = {E_-E + )/2. 



(2) 



This relation may be understood as follows. Consider 
an electron which hops between two atoms of the same 
atomic energy with hopping integral t. Then, the energy 
of the even parity state is t and the energy of the odd 
parity state is —t, the energy difference between the two 
states is (E + — = 2t. In the ion problem, the 
singlet spin gives an additional sign change, as does Eq. 
(§• 

The H 2 ion is a three-body problem and therefore no 
exact analytic solution is available. Yet a suitable ap- 
proximation can still give a reasonably good ground state 
energy. We construct the lowest energy states of different 
parities using the single site states of H~ ion, Eq. (Q) , and 
the neutral H-atom. Let 4>{a) be the hydrogen ground 
state wavefunction (the Bohr atom solution) of the ath 
electron (Fig. ||). Then the three electron states of the 
with odd(— ) and even(+) parities are given by: 



where 



*± = <!>,.., ± '!',.,. 



$^•=4^(1,2)^(3)]. 
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A is the antisymmetric operator to assure the anti- 
symmetry of the wavefunction when two electrons 
are interchanged. The corresponding energies of 'f ± are 
given by 
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where h is the Hamiltonian for the system, 
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h = Y(^- — -—) + T — 
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In Eq. (H), /2m is the kinetic energy, r a and r' a are 
the distances of the ath electron to the two protons re- 
spectively, and r a/ 3 = \r a — rg|. It is convenient to split 
h into two parts, h = ho + h', with ho the Hamiltonian 
without the interaction between the sites i and j , 
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ft-o approaches to h at the limit of large inter-proton dis- 
tance, h' describes the interaction energy between the 
H~-ion and the H atom. 

Substituting 4>j, ,- into Eq. (||) in terms of fa and fa , we 
find the two numerator terms to be the following combi- 
nations of integrals of ho and h' , 

(<Z> iJ \h\$ i j) = 3(Eo + E 1 )-3(Eoa2 + b 2 ), (8) 
($j,i\h\$ij) = S(E Q a 3 + b 3 ) - 3(E oai + 6i), (9) 

where E' s and a' s are defined as, 



very successful on capturing the most essential physics 
in a single H^~ ion. This also justifies our estimate for 
the electron hopping integrals in the H~ ion lattice. The 
present calculations give a higher ground state energy 
for d < 4 a.u. than that of Taylor and Harris. This 
is understandable. For smaller distances d, our method 
constructing a wavefunction based on individual H and 
H orbitals becomes poor. However, the distance between 
hydrogen ions in lanthanum hydrides is always greater 
than 4 a.u., the region of our interest. Therefore we are 
confident in our estimate of the hopping integral. 

Using Eqs. (||) and (|ll"|), we find that the hydrogen 
hopping integral t can be written in terms of the integrals 
in Eq. (foh as, 



E Q = ((12) i 3 J |h |(12) i 3 i ) 
a x = ((12) i 3 J |l i (23) i ), 
a 2 = ((12)^1(23)4,), 
a 3 = ((12)^13,(12),), 



E, = ((12) i 3 i |/ l '|(12) i 3 i ), 
h = ((12) i 3 j |/i / |l i (23) j ), 
b 2 = ((12) i 3 j |/i'|(23) i l J -), 
6 3 = ((12) i 3 J -|/ l / |3 i (12) J -). (10) 



For simplicity, we define | (12)^3-,) = |^(1, 2)0,(3)). Note 
that Eq is simply the sum of the energies of the indepen- 
dent H _ ion and H-atom. Eq. (0) now becomes, 
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Ei-b 2 ± (63 - h) 
(1 - a 2 ) ± (a 3 - 01) ' 



(11) 



These integrals are calculated numerically for a few val- 
ues of the inter-proton distance d and are listed in Table Q. 

For H^ ion, the ground state is of odd parity. The 
energy of H^~ ion we obtained in the present approach 
compares favorably with the best estimates reported pre- 
viously using different methods In Fig. [|, we 
plot our values of E_ as a function of inter-proton dis- 
tance d. Fischer-Hjalmars |25j used a similar approach, 
but ignored the correlation term in the Chandrasakhar 
wavefunction probably due to the limited computational 
power in the 60's. Our energies are much lower than 
theirs. This indicates that the electron correlation has 
a significant contribution to the ground state energy of 
the H^ ion. It comes as a surprise that our result is even 
somehow better than that of Taylor and Harris for d > 4 
a.u.. Taylor and Harris (24| used a rather complicated 
wavefunction which involves many linear and nonlinear 
variational parameters. Their parameters are distance 
dependent and variational procedures were carried out 
for each of d. The comparison with their result indicates 
that the Chandrasakhar wavefunction that we used is 



, Cgl " b 2 )(ai - a 3 ) - (61 - b 3 )(l - a 2 ) 
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giving values for the hopping integral between the neigh- 
boring H tet and H oct (see Fig. 0), t 2 = -0.748 eV 
(d — 4.58 a.u.), and between the two nearest H tot atoms, 
ti = -0.523 eV (d = 5.29 a.u.). 

In the above calculations, we focus on the H^~ ion, and 
have neglected the Madelung potential of other H~ ions 
and metal ions in the lattice. These ions create an electric 
field on the H^ ion under consideration, which reduces 
the amplitudes of hopping integrals. This effect is often 
called a crystal field effect, and will be discussed next. 



C. Crystal field effect 

In the above estimation for the hopping integrals, we 
have neglected the La 3+ ions as well as surrounding H~ 
ions. In an ionic picture, the La 3+ and H _ ions generate 
a crystal field at each H site. Therefore the Hamiltonian 
in Eq. @ should be modified to include the Coulomb in- 
teractions between the electrons in the H^~ ion under con- 
sideration and all other H~ ions and La 3+ . We will treat 
the crystal field effect as a perturbation, whose Hamilto- 
nian is given by 



Z(R)e 2 



R 



r 2 



Rl' 



(13) 



where the sum runs over the La 3+ , and all other H 
ions. Z(R)e represents the charge of an ion located at 
R ( Z(R) = 3 for La 3 + ion, and Z(R) = -1 for H" 
ion). Then the first order correction to the energies E± 
is given by 



5E± 



E[ T b[ 
lT«i ; 



with 



^=((12)^1^1(12)^), V L = ((12) i 3 j \h" 



(14) 

1,(23),). 
(15) 
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Accordingly, the change of the hopping integral, St 
{5E- - SE+) /2 is found to be 



St = - 



E[ai - b[ 



(16) 



In the calculation of St, we shall use a summation tech- 
nique similar to the usual calculation of the Madclung 
constant to keep the charge neutrality of the summed 
ions. For a reason to be explained later, we consider 5t\, 
the change of the hopping integral between the nearest 
neighbor (n.n.) tetrahedral hydrogens. We consider the 
midpoint of the H^~ in question as the origin (Fig. (0)), 
and calculate St\ from the contributions of all the other 
ions within a cubic cell centered at the origin. As the size 
of the cubic cell increases, 8t\ changes monotonically and 
saturate. In Table ||, we list 5t\ as the size of the cubic 
cell increases. 8t\ starts to converge for a 4 x 4 x 4oq cu- 
bic cell and reaches a value 8t\ = 0.230 eV at 6 x 6 x 6<2q 
(See Table @). 

Note that t\ is negative , and 8t\ is positive. The crys- 
tal field reduces the magnitude of the hopping integral. 
This is consistent with our intuition. The net effect of 
surrounding ions is from positive charged ions, because 
of the charge neutrality (H^~ has a net charge of — e). 
The positive charged ions (La 3+ ) attract the outer elec- 
tron of a H~ ion, and reduce its hopping amplitude to a 
neighboring hydrogen. 

We thus estimate the hopping integral in the presence 
of the crystal field, to be t\ — —0.293 eV. This is about 
60% of the value without the crystal field. We see that 
the crystal field significantly reduces the hopping inte- 
gral. However, the estimate for t2 is more complicated 
because of the asymmetry of the crystal with respect to 
the tetrahedral and octahedral hydrogen atoms. Assum- 
ing the same percentage reduction for ti, we estimate 
?2 = —0.419 eV. Note that our estimate is based on an 
ideal ionic approach within which screening effects are 
neglected, so that the actual reduction of the i's is ex- 
pected to smaller. 



III. MICROSCOPIC MODEL FOR METAL 
HYDRIDES 

A. Model Hamiltonian 

Having a better understanding of the electron correla- 
tion in H _ and its effect on the electron hopping matrix 
between the H ions, we are ready to proceed to an ap- 
propriate microscopic model for the rare earth hydrides. 
We introduce a large [/-limit Anderson lattice Hamilto- 
nian to model the system. In the tight binding repre- 
sentation, the Hamiltonian, in the second quantization 
language, is given by: 



H — Hh + Hlc, + H„ 



(17) 



where the three terms represent the hydrogen, the lan- 
thanum, and their hybridization respectively. The hy- 
drogen part is 



h - (h f h 



h.c. 



(18) 



where i sums over all the occupied H atoms, and 
neighboring pairs. hi tS is the destruction operator for an 
electron of spin s on the H site i. There is a constraint 
for electrons on each H site i, 



E 



h] Mi , > 1. 



(19) 



Eq.(|l9|) is to exclude the empty electron state at any 
hydrogen site. We will discuss this point further below. 
In Eq. (|T8]) , is the atomic energy of the outer electron at 
site i, which is —0.7 eV in free space and will be modified 
in crystal, e^ 1 = e t at the tetrahedral site, and — e at 
the octahedral site. & are the hopp ing i ntegrals between 
two H sites as estimated in section [IB . The lanthanum 
part of the Hamiltonian is given by 

///.,• = E <':"< f ,. J ■ + E e d ( d L, A/m + 



(hi) 



(20) 



where di tCX ,s destroys an electron of orbital a and spin 
s on the La-site i. We shall only include the five La-5d 
orbitals labeled by a. They are closest to the chemical 
potential and therefore most relevant. 

The La-6s orbital will be neglected since its energy 
levels are high above the fermi energy. e^ a denotes the 
atomic energy of orbital a at the the La-site, tf'? is the 
hopping integral between orbital a at site i and orbital 
(3 at site j. 

H m ix describes the electron hybridization between H 
and La sites, 



Hrr 



E 



v,. 



(<«,- 



h.c 



(21) 



where Vi. a ,j is the hopping integral between La site i of 
orbital a and H-site j. 

We now discuss on the physical meaning of the con- 
straint Eq. (|l9|) . As we examined in the previous sec- 
tion, the two electrons in the H _ ion are very different in 
energy. The outer electron has a binding energy 0.7 eV, 
while the binding energy for the inner one is 13.6 eV in 
free space. At low energies, the inner electron is always 
occupied, and only the outer electron is mobile. The con- 
straint (|l^) is a mathematical description of this physics. 
If we define a H~ ion as a vacuum, namely a fully filled 
Is shell, then a neutral H atom is a single-hole state, 
and H+ is a double-hole state. The constraint Eq. ( pj| ) 
prohibits a doubly occupied hole configuration at any H 
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site. Our model ( |l7| ) is a large [/-limit Anderson lattice 
model, where U is the Coulomb repulsion between two 
holes on the same H site. 

The important difference between inner and outer elec- 
trons in hydrides was also discussed by Eder et al. 
These authors describe hydrogen as a breathing atom, 
whose radius is much larger for H~ than for H. Eder et 
al. proposed an Anderson lattice model, where the hop- 
ping integrals depend on electron occupation number on 
the H site. They carried out an impurity-like calculation 
to examine the stability of the Kondo-like singlet con- 
figuration similar to that proposed in our earlier work 

We now apply Eq. ( |17| ) to LaH^. In LaH^ the lan- 
thanum sites are always fully occupied while the occu- 
pation of hydro gen sites depends on x. The properties 
of Hamiltonian (jl7l)-([jl"|) also depends on the parameter 
values. The effective hopping integrals t^j b etween two 
H sites have been estimated in section II B. The other 



and Vi ta j can be represented by a few bonding integrals. 
They are given by 



parameters may be extracted from the local density ap- 
proximation calculations, which will be described in the 
next section. 

The Hamiltonian (|l7]) with the local constraint is a 
many-body problem. We will solve the model for differ- 
ent values of the H concentration x using the Gutzwiller 
method |^8|,^9| , which has been well established in study 
of heavy electron and mixed valence compounds. 



B. Parameters fitting from LDA calculations 

The LDA does not treat the electron correlation prop- 
erly. However we expect the LDA to give reasonably good 
estimates for other parameters in the hydrides. We car- 
ried out LDA calculations for LaH3, used a tight binding 
model to fit to the LDA results and extracted parameters 
for model (fL7|). For simplicity, we consider only hopping 
integrals between the nearest neighbor (n.n.) tetrahedral 
H-sites ^2), between the neighboring H to t and H oct sites 
(ti), between the n.n. H and La sites, and between the 
n.n. and the next n.n. La sites. Since the lattice con- 
stant changes insignificantly for 3 < x < 2 in LaH^, we 
shall assume the parameters remain unchanged. 

In Fig. [7] we show our LDA results for the band struc- 
ture of LaH3 and LaH2. Other LDA calculations gave 
similar results 
identical to Eq 




We use a tight binding model, 
but with no constraint Eq. to 
fit Fig. ^(a) for LaH3. The fitted band structure is plot- 
ted in Fig . ||(a), and the fitting parameters are listed in 
Table |H4 In the fitting, we pay most attention to those 
states near the chemical potential. The same set ofpa- 
rameters fit the LDA bands of LaH2 as well, see Fig. 0(b) 
and Fig. |(b). 

Since all s-s, s-d and d-d orbital bondings can be ex- 
pressed as a linear combinations of cr-bonds (angular mo- 
mentum along the bonding axis m = 0), 7r-bonds (m = 1) 



j-h xrh — h 

i.j ssa 1 



ttf=a a ^l,m,n)V s L d l- h , 
Vi^ a .j = b£\l,m, n)V d da + b£\l,m,n)Vddir- 



(22) 



where V^ s ~ h {V^ J ~ h ) is er-bonding integrals between H 
Is and H Is (La 5d) orbitals, while Vdda and Vm-k & re 
cr- and 7r-bonding integrals of La 5d to 5d orbitals, re- 
spectively. The contribution from <5-bonding is generally 
small and is ignored. The coefficients a a ^ and 6 Q 's are 
functions of the direction cosines (^,m,n) from site i to 
site j. These coefficients have been derived by Slater and 
Koster |30| based on the symmetry. The number of fit- 
ting parameters is now reduced since there are only a few 
bonding integrals. It is worth to note that the H hop- 
ping integrals obtained from LDA, V"^ = —0.79 and 
Vg~J — —0.4, are larger than that calculated from the 
local model, t 2 = —.419 and t\ = —0.293. It justifies 
that the electron correlation plays an important role in 
the H hopping integrals. Therefore, we will use t\ and t 2 
for the H hopping integrals in the future calculations. 

Although there are some discrepancies with the LDA 
band structure for high energy La bands, the main con- 
cern here should be the position of the lowest La bands 
and the highest H bands which are responsible for both 
the transport and optical properties. Bear in mind that 
our motivation is to reproduce the essential qualitative 
features rather than a very accurate description and 
hence the fitting may not be unique. We have carried 
out the same fitting procedures for Dekker's LDA band 
calculation jljj. Even though it generates a different set 
of parameters, the general features like the opening of 
band gap when introducing electron correlation and the 
localization of the vacancy and impurity state, which will 
be discussed later, are consistent with the current set of 
parameters. Therefore, we are confident in the use of this 
parameter set for our future calculations. 



C. Gutzwiller method for the many-body 
Hamiltonian 

Following Rice and Ueda we shall use the 

Gutzwiller method to study the many-body Hamiltonian 
(|i"7j). The method is a variational approach. The con- 
straint condition of Eq. (|l9|), to exclude state at any H 
site, is replaced by a set of renormalization factors. 

We consider a variational ground state wavefunction 
|\&) which is of the form 



I*) = P\^o 



(23) 



and 5-bonds 



2), the hopping integrals t 



where |\& ) is a state of the tight binding Hamiltonian 
with no constraint, and P is the projection operator to 
project all empty electron states on each H site, 



G 



P : 



n 



(24) 



where i runs over all H sites and 



The ground state energy is given by 

(9\H\V) 



E„ 



(H) 



(25) 



To calculate (H), we use the Gutzwiller method to com- 
pute the average values of each term in H in the state 
\ty). The Gutzwiller method is a static approximation, 
and the average of a local operator Q in the state \ty) 
is related to its corresponding average in the |\l/ ) by a 
renormalizcd numerical factor gQ, namely, 



(Q) =9q(Q)o, 



(26) 



where (Q) and (Q) are the average in the states |*) and 
\^ ) respectively. 
In the present case, we have 



i,A>) = V9i9j( h ls h i,s)c 



(h 



(27) 



The Gutzwiller factor g is calculated by counting the 
number of possible configuration in \^f) and in |W D ), re- 
spectively. The value of is 



1 



-, hole 



1 -n 



holi 



/2' 



(28) 



where 



-,/io/c 



is the occupation number of holes at the H 



site i. Here we define a hole as an electron vacancy in a 
H - ion and therefore 



hole 



The constraint Eq. (|T^) now disallows doubly occupied 
hole configurations on H sites. This leads to an effective 
Hamiltonian for the lanthanum hydrides, 



H eff - H h + H La + H'raix ~ M ^ h i,s h i 



(29) 



In the above equation, Hi, a is given by Eq. pOj ) , and Hi 
and H' ix have the same form of Hh and H m i X in Eq. Jig ) 
and The constraint condition is released, and the 

hopping integrals are renormalized to values 



(30) 



H hopping integrals ij ■ are the values obtained from our 
local model including the crystal field effect. Further- 
more, the atomic energies e' 1 are shifted to — ej,, where 
£b is the binding energy of the outer electron in a H~ 




ion. We use the free space value e;> = 0.7eV. Note that 
the outer electron in LDA is unbound (2?J . /u Eq. (H|) is 
an effective "chemical potential" J3lj| . The value of /i is 
determined by minimizing the total energy of the system 



Eg = (H eff +nJ2 h Lh t 



(31) 



In the lanthanum hydrides, the total energy is found to 
be insensitive to values of \i ~ 0, so we set \i = 0. This 
is primarily due to the fact that e h is lower than e a . 

To compare with heavy fermion systems, we may per- 
form an electron-hole mapping for the lanthanum hydride 
model Hamiltonian, e h — * — e , e a — > — e a , and the con- 
straint is to project out doubly occupied hole states on 
H sites. In the hole notation, the "f level in lanthanum 
hydrides is above the conduction bands. Therefore the 
hydrides are in the mixed valence, instead of the Kondo 
limit of the Anderson lattice model. 

H e ff is a single particle Hamiltonian, and can be 
readily solved by analytic or numerical means. The 
Gutzwiller factor gi, given by Eq. (p8|) , satisfies a self- 
consistent equation for the electron occupation at H site 



n't = (n?) = (nf) . (32) 

Since < g < 1, the hopping integrals are then renor- 
malized to smaller values. Hence the hydrogen band 
width is expected to be further reduced as a consequence 
of the correlation. In general g is position dependent. In 
the periodic cases such as in LaH2 and LaH3, g depends 
only on whether it is on the tetrahedral or octahedral site. 
For the inhomogeneous impurity states, g are strongly 
site dependent. 

Note that this occupation related correlation is dif- 
ferent from the intr asite electron-electron correlation we 
discussed in section II B. There the interaction between 
electrons on the same H~ ion was considered, and the re- 



sult is the correlated wavefunction of two electrons on the 
same site. Here in this section we consider the correla- 
tion effect on the intersite hopping, which is a constraint 
on the every H-site. 



IV. ELECTRONIC STRUCTURE OF LAH 2 AND 
LAH 3 

We are now in position to calculate the electron spec- 
trum of lanthanum hydrides. In this section, we consider 
LaH2 and LaHj. 

In the case of LaH2, all tetrahedral hydrogen sites are 
occupied, and all octahedral H sites are unoccupied. Be- 
cause of the equivalence between the two tetrahedral H 
atoms in the unit cell, we need to introduce only one 
Gutzwiller factor, g t . The energy spectrum, or the band 
structure, obtained from the H e ff, is plotted in Fig. ||. 
The two lower bands are primarily H _ bands, and are 



7 



fully filled, leaving one electron per unit cell in the bands 
of primarily La-5d characters. The Fermi energy is within 
the La-5d conduction bands, so that LaH2 is a metal as 
expected. The correlations in LaH2 do not change this 
qualitative property. 

In LaH3, all octahedral H sites are also occupied. We 
use two independent Gutzwiller factors gt and g corre- 
sponding to the tetrahedral and octahedral sites respec- 
tively. The Gutzwiller factors are found to be g t = 0.78, 
g = 0.70, which indicates the reduction of the hopping 
parameters between two H sites and between H and La 
sites. The electronic structure for LaH3 is plotted in 
Fig. 0. 

We find that LaH 3 is an insulator. The optical gap 
is found to be 1.5 eV or 2.1 eV depending on whether 
the crystal field is included in the estimate of the H-H 
hopping integrals as discussed in section II. The opening 
of the energy gap at the chemical potential is primarily 
due to the strong electron correlations in the H~ ions, 
which reduces the H~ band width. The large electron- 
electron Coulomb repulsion on the H ion also restricts 
the electronic state to exclude the empty electron state 
at any H-site, which further reduces the H - band width. 
The energy of the conduction band at X-point seems too 
low, and is presumably due to the inaccuracy in the tight- 
binding fit, see Fig. 0(a) and||(a). 

The present method is an improvement of our previ- 
ous work jl6j. In that paper, we focused on the H bands 
and neglected the La bands. That is to assume that all 
the three lower bands, which is of primarily H-ls char- 
acter and mixed with the La-5<i orbitals, are purely H 
bands. Mathematically, the previous model is equivalent 
to the effective Hamiltonian Eq. (|2^) but without the 
renormalization factors (i.e. gi = 1). The more sophisti- 
cated method presented here predicts a larger gap, closer 
to the observed value. 

It would be interesting to study the band gap in LaHs 
as a function of the lattice constant a . As the lattice 
constant decreases, the electron hopping matrix elements 
between H-ions increase, and thus the H - band width 
increases, possibly leading to the reduction of the gap 
and eventually to the closing of the gap. Therefore we 
expect a transition from insulator to a metal in LaHs 
under pressure. With the knowledge of the distance de- 
pendence of the hydrogen couplings t\ and £2, crystal 
field effect (cx 1/ao) and the orbital coupling parame- 
ters (VddTv and Vdda oc V sda <x 1/a^/ 2 ) fhj, we can 
calculate the evolution of the band gap as a function of 
a . We find that the band gap closes and LaH 3 becomes 
metallic when oq decreases by 11% to 0.89a o . (Fig. 11). 



V. IMPURITY AND VACANCY STATES 

In the previous two sections, we argued that electron 
correlation is the missing element from the LDA calcula- 
tions, which failed to predict the insulating properties of 



the trihydrides. By introducing an electron correlation 
on ions, we have found that LaH3 is an insulator, 
consistent with experiment. However, more interesting 
physics lies in the mechanism of the metal-insulator tran- 
sition in LaHa; . We can understand the transition by ap- 
proaching it from either metallic (x — 2) and insulating 
(x = 3) end points. We investigate the former first. 



A. An extra H atom in LaH2 

The introduction of a neutral H atom into a H oct site 
in LaH2 creates a s=l/2 magnetic impurity, which cou- 
ples to the conduction electron spins. Let us study this 
impurity problem more carefully. First we consider the 
atomic limit and neglect the hybridization between the 
H-ls and La-5d orbitals. Then a conduction electron 
(La-5d) moves to the neutral H atom at H oct site to form 
a H~ ion. Including the hybridization between H and 
La orbitals, the outer electron in the H~ ion will fluctu- 
ate between the H oct site and the surrounding La sites. 
Since H _ ion is a spin singlet of two electrons, the cou- 
pling between H oct and La orbital is antiferromagnetic. 
The effective Hamiltonian for the extra H atoms at H oct 
is an antiferromagnetic Kondo model, 



(33) 



k,a,s 



where i runs over all the occupied H oct , S and s are the 
electron spins of the neutral H atom and of the conduc- 
tion electron (annihilation operator dk,a,s with energy 
e Q (k)) states. Because of the symmetry, only the two 
e g La-5d orbitals are coupled to the ls-H oct electron di- 
rectly. The exchange coupling will be large, J ~ 0.7 eV is 
in free space To estimate J in the hydrides, we use 
a local model including six neighboring La atoms around 
the H oct site (Fig. [j^). The dynamics is given by a 2 x 2 
matrix Hamiltonian 



tact V6V\ 



(34) 



where e oc t is the atomic energy of the 1-s H oct and V = 
VJdo~° 1S t ne hybridization between the H oct and each La. 
id is the energy of the combination of the six 5d-La-e ff 
orbitals around the H oct , denoted by 



\ im p) = ~7g X! 
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(35) 



at the La 



where |e* ) is the e g orbitals as shown in Fig. 
site i. Using the Hamiltonian H^ a in Eq. (20), we find 
the energy of \imp) is to be given by 

e d = (imp \H La \ imp), (36) 

= e'd + €k in ■ 

tKin is the kinetic energy of the electron in the \imp), 
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Vdda 



iVddll, 



(37) 



where t eg ~ ea is the hopping integral between two La-e g 
states whose orbitals point to the central H oct . From 
Eq. (S3), we find the spin-singlet state energy to be 



£s=0 



(eoct + 



e d 



2W 2 



(38) 



Using the parameters fitted to the LDA, e oct = —3.3 eV, 
e d = -0.7 eV, V = -1.2 eV, we have e s=0 = -5.21 eV. 
This energy should be compared with the chemical po- 
tential (fi = —0.62 eV, see Fig. ^) for a spin triplet state, 
in which an electron with the same spin is placed at the 
Fermi level conduction band. Therefore, we arrive at a 
value 



J = (jl- e s=0 = 4.59 eV. 



(39) 



In this Kondo description conduction electrons are cap- 
tured by the neutral H atoms at H oct sites to form tightly 
bound singlets. In this description the material with all 
octahedral sites occupied, LaH3, is viewed as a Kondo 



insulator with a large band gap. 
is consistent with Eder et al. [11 
the insulator is expected to be reduced from the value in 
Eq. (p9|) obtained within an impurity calculation. 



Our conclusion here 
The actual gap in 



B. Localized vacancy state in LaH3 

Next we examine the transition starting from insulat- 
ing LaHV The removal of neutral H atoms introduces 
vacancies at octahedral sites H^, t which donate an elec- 
tron to the conduction band. In a conventional semicon- 
ductor such as in the phosphorus-doped silicon (Si:P), 
the impurity state is described by an effective mass the- 
ory, and the result is a hydrogen-like bound state with a 
large effective Bohr radius a* B of order ~ 100 A due to 
the light effective electron mass and the large dielectric 
constant. The critical impurity concentration S c at which 
the system becomes metallic is given by Mott criterion 
5 C ~ (1/a.g) 3 , where I is the inter-atomic distance. Be- 
cause of a* B >• I, S c is very small (10~ 3 for Si:P). The 
vacancy state in L&H3-S is very different, however. Ex- 
perimentally, it is found that the semiconducting states 
extend to S c = 0.25 for LaH3_,j, and to even larger value 
for yttrium hydrides |i|,[Tl| . This suggests a very localized 
impurity state. In our previous short paper we proposed 
the vacancy state in LaH3_,5 to be a very localized La- 
5d-e g orbitals centered at the vacancy with s-wave sym- 
metry, and we used perturbation theory to estimate its 
energy in support of our proposal. Here we calculate the 
lowest vacancy state and further support our previous 
vacancy model. 

Let us begin with a qualitative argument and a sim- 
ple estimate. The La-5d-e ff and H oct electrons hybridize 



strongly to form bonding (mainly H) and antibonding 
(mainly e g ) states. The underlying physics may be sim- 
ply described by the local model of a single H oct atom sur- 
rounded by six La atoms (see last section). The bonding 
state energy is given by Eq. (38), and the anti-bonding 
state energy is 



e a b 



(eoct + id) + \J {e ct -idf + 24V 2 



(40) 



A vacancy of H oct breaks the bonds locally and the e g 
electron becomes locally non-bonding, which has much 
lower energy than the antibonding e g states away from 
the vacancy. Therefore, in addition to serving as a posi- 
tive charge center as in conventional n-type semiconduc- 
tors, H^, t creates a potential well for the e g state electron. 
The latter is non-perturbative, and is responsible for the 
unusual concentration dependence of the semiconductor. 

Consider a s-symmetric octahedral e g state (5-state 
hereafter) around a H y t as shown in Fig. [l| The H v t 
vacancy breaks the bonds on the octahedron and reduces 
the antibonding energy of the 5-state. This 5-state is 
very localized because the neighboring octahedral s-state 
is antibonding and has a much higher energy. We can 
estimate the depth of the potential well, V$. The vacancy 
state (the non-bonding state) energy is 



Evac — fzd- 

Therefore V — E a b — E n b is given by 



e oct -e d +y {e ct - id) + 24V 2 



(41) 



(42) 



The numerical value of the potential well is V = 1.9 eV. 

From Eq. (fll|), we estimate E n i = —0.7 eV, below the 
lowest La-5d-t2o band which has energy zero as we can 
see from Fig. |lCl|(a). 

In the above estimate for the energy of the localized 
vacancy state, we have not included a Coulomb attrac- 
tion between the H oct vacancy and the extra electron. 
This attraction will further lower the vacancy state en- 
ergy. We conclude that the H oct -vacancy electron state 
is highly localized and it is well within the band gap. 
The localized vacancy state here is similar to the carbon 
vacancy state in TiC, where the C-vacancy leads to a 
well localized impurity states p3] |. We emphasize that 
the potential well generated by the H oct vacancy is pri- 
marily short-ranged, different from the ideal long-range 
Coulomb force. The effective radius of the impurity state 
is only half of the lattice constant. 

The existence of localized vacancy state can be ver- 
ified by a finite size calculation for a single H oc t va- 
cancy in LaH 3 , an improvement of the local model dis- 
cussed above. The system is described by Hamiltonian 
(|l7|)-(p8|) with all but one H oct o ccupie d. We use the 
Gutzwillcr method given in section [II C to study the va- 
cancy state. The vacancy breaks the translational invari- 
ance; the Gutzwiller factor gi now depends on the site, 



9 



and has only crystal group symmetries with respect to 
the vacancy state. Because of its localized nature, we can 
study the vacancy state in a finite-size system. The prob- 
lem then can be solved explicitly using the Gutzwillcr 
approximation. We consider a small cubic of 3 x 3 x 3a^ 
lattice in the calculations. We choose the vacancy site 
to be at the center of the lattice. An S'-state around the 
vacancy site is found, as we expected. All the surround- 
ing La atoms have their 5c? e g orbitals pointing towards 
the vacancy and have the same weighting factor. The tig 
orbitals, however, have no contribution in the S state. It 
is found that the S state has an energy E vac = —0.38 eV 
below the 5d La conduction band at T point which is a 
tig state and is decoupled from the vacancy state. There- 
fore, electrons will occupy the less energetic vacancy state 
instead of the conduction band, implying an insulating 
ground state for x < 3. When more vacancies are added, 
more localized states are formed within the band gap and 
the material will become metallic when localized states 
start to substantially overlap with each other. 

In Table [v| we list the probability of the vacancy state 
at different atomic sites in the lattice. As we can see, 
the state is largely distributed in the six neighboring La 
sites with the probability of 58%. This clearly shows its 
localized nature. 

The localized nature of the vacancy state in LaH^ is 
consistent with the temperature dependence of the d.c. 
resistivity data at room temperatures, which has a tem- 
perature depedence consistent with variable range hop- 
ping . Localized states are a prerequisite for variable 
range hopping. In the next section we will discuss the 
electronic band width of the vacancy states. 



C. Effective electron hopping integrals between two 
vacancies 

In LaH3, octahedral hydrogens are not fixed in posi- 
tion but able to move around by diffusion. Electrons 
in the vacancy state, of course, are strongly bound by 
the Kondo like effect as discussed before. Therefore, a 
calculation of the vacancy-vacancy interaction and the 
effective electron hoppings between the two H^ ct is nec- 
essary to the electronic band of the vacancy states. We 
estimate here the effective hopping integrals between two 
Hoct sites. 

The effective hopping integral between the two 
H oct vacancy sites is given by 



t v . v = |i?'|*2>, 

where H' is the inter-atomic Hamiltonian given by H in 
Eq. (|l7|). I^i) and are the electron state of the two 
vacancy sites. 

Note that the vacancy states \*f? a ) are basically linear 
combinations of the nearest 6 L&-5d-e g orbitals (Fig. 14). 
For simplicity, we model the vacancy state by the local 



description as shown in Fig. |l3|, and only the n.n. cou- 
pling between La-sites in H', and hence Vdda and Vdd-K, 
will be considered while contributions from further dis- 
tant orbitals is negligible. We apply the fitted values of 
Vdda and Vddir from section [II. Table VI lists the val- 



ues of hopping integral of vacancy states separated by 
different distances. 

It is interesting to notice that the value of the hopping 
integral oscillates when the distance increases. A more 
careful investigation tells us that the hopping of the va- 
cancy state within the same cubic sublattice is larger in 
value and negative in sign. In La hydrides H oct sites 
forms a face-center cubic lattice, which can also be di- 
vided into 4 interlocking simple cubic sublattices. Our 
finding shows the hopping within the same sublattice is 
more probable than between different sublattices. 

The precise nature of the transition in optical proper- 
ties is not yet well understood. Since the vacancy states 
are highly localized, the vacancy concentration at the 
metal insulator transition of a random vacancy distribu- 
tion is quite high. Hence the mobile carrier density is 
high. This explains why the metallic phase is a good 
metal with high reflectivity. 

One aspect of the hydride system is unique, namely 
the existence of an order-disorder transition among the 
vacancies which enables one to separate the effects of dis- 
order and Coulomb interaction in the electronic band of 
vacancy states. Usually donors and acceptors are fixed in 
doped semiconductors so that this is impossible. The dis- 
ordered vacancy phase is insulating with variable range 
hopping among the localized vacancy states. Interest- 
ingly Shinar et al. jllj| reported a transition to an or- 
dered vacancy phase at temperatures around 250°K for 
concentrations x ~ 2.8 — 2.9 and metallic d.c. conduc- 
tivity in the ordered phase. This shows that Anderson 
localization is responsible for the insulating character of 
the disordered phase. However in the ordered phase there 
must be a Mott transition to an insulator as the vacancy 
density is lowered. It would be very interesting to ex- 
amine this region and look for the Mott transition in an 
ordered vacancy sublattice. 



VI. PHYSICAL PROPERTIES OF INSULATING 

LAH 3 _i 

In order to study the electronic structure of metal hy- 
drides, various measurements, such as photoelectron and 
optical spectroscopy, have been performed for metal hy- 
drides with different hydrogen concentration. From these 
experiments [j36| interband excitations can be measured 
and compared to theoretical band structure. Their data 
can be used to verify the accuracy of theoretical predic- 
tion for the electronic structure. Most theoretical investi- 
gations have been done on the dihydrides and trihydrides, 
while hydrides with intermediate hydrogen concentration 
are difficult to study. In this section we will discuss some 
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physical properties predicted from our theory in connec- 
tion with the existing or further experiments. 



A. Magnetic properties of LaH 



3-<5 



The vacancy state i s a localized object of spin-1/2 



as discussed in section 
Lari3_,5 is an insulator for small <5 



VB . Near the trihydride phase, 
The vacancies have 



a tendency to be singly occupied because two electrons 
on the same vacancy site repel each other. The magnetic 
susceptibility is expected to follow a Curie law with a 
linear inverse temperature dependence. 



B. Density of states 

Peterman et al. |36| studied the composition- 
dependent electronic structure of LaH^, 1.9 < x < 2.9, 
using photoelectron spectroscopy with synchrotron radi- 
ation (10 eV < hv < 50 eV). Fujimori et al. |}7| on the 
other hand studied the electronic structure of yttrium 
hydride by x-ray photoemission spectroscopy. Here we 
compare our calculated densities of states with the ex- 
periment on lathanum hydrides, as shown in Fig. [l5| 

We performed a finite size calculation with lattice size 
of 15 x 15 x 15ciq based on the model Hamiltonian (17). 
Numbers of states in the band structure of LaBi3 and 
LaH2 are counted as a function of energy. The dis- 
crete densities of states were smoothed out by replacing 
delta peaks with Lorentzian distribution curve of width 
r = 0.3 eV. Our results for LaEb and LaBi3 are sim- 
ilar to the pre vious theoretical calculations by Gupta 
and Burger p6| , p^8| . In Fig. the pronounced two-peak 
structure is associated with the flat region of hydrogen 
bands near symmetry points X, L, K (refer to the band 
structure Fig. ^| and [n]) while the small bump at the 
zero energy in LaH2 is associated with the La band near 
symmetry points W and K. Experiments showed that the 
small bump shrinks as hydrogen concentration increases 
and eventually disappears when the sample approaches 
the trihydride. However Gupta and Burger expected a 
metallic state for LaH3 from their LDA calculation which 
ignored electron correlation. In our calculations, the den- 
sity of states goes to zero at —2.2 eV, indicating an 
insulating behavior for LaH3. However, the calculated 
widths for LaH2 and LaH3 are considerably smaller than 
the experiment results, a discrepancy similar to that in 
the LDA . The experimental data also shows a shift of 
the lowest energy peak when the concentration increases. 
The positions of calculated peaks do not significantly de- 
pend on the concentration and does not agree to the ex- 
periment. The origin of these discrepancies is not clear. 
Since our theory uses some parameters extracted from 
LDA, similar features except for the band gap may be 
expected from both theories. 



C. Optical Conductivity 

The optical absorptivity spectra for samples of LaH 2 . 
and NdH-j were also measured by Peterman et al. |3(| 
at 4.2K at near-normal incidence. From the absorption 
data, they deduced the real and imaginary part of the 
dielectric constant, and hence the optical conductivity 
by using Kramers-Kronig analysis. They found a rela- 
tively broad feature for LaH2.87- However, their samples 
were polycrystals and the use of Kramers-Kronig analy- 
sis might enlarge the data uncertainties. On the other 
hand, a more recent and accurate measurement reported 
by Griessen et al. |59| sheds new light on the optical con- 
ductivity curve. The optical transmission spectra of the 
insulating phase YH3_^ were measured as functions of 
the photon frequency Tito and of hydrogen vacancy con- 
centration 8. The effect of the vacancy appears to reduce 
the overall transmission spectra quite evenly between 
huj = 0.5 eV and 2 eV. In our theory, the (5-dependent 
conductivity <j{lo) at hu < 2.8 eV mainly determined by 
the optical transition from the vacancy state to the con- 
duction bounds. The transition energy from the valence 
bands to the vacancy state is larger because of the larger 
energy difference between the two states and because of 
the Coulomb repulsion of the doubly occupied electron 
states on the same vacancy site. Since the vacancy state 
is highly localized, the transition matrix largely depends 
on the density of states of the conduction bands, which 
is expected to lack pronounced feature. 

Within the dipole approximation, the real part of con- 
ductivity contributed from the vacancies is given by 

er (ui) — const x uj | (n \x\ vac) | 2 5 (E n — E vac — Tilj) 

n 

(43) 

where \vac) is the vacancy state. The sum is over all 
conduction states In). We shall use a local model to 



describe \vac), as shown in Fig. 13 



— ^ a Q (R_L a )|RL Q , do,), 



(44) 



where R_L a denotes the position of the 6 La atoms and 
d a the 5d-e g orbitals, and a a (Ri n ) are the wavcfunc- 
tion amplitudes. The conduction bands are mainly 5d- 
La character mixed with ls-H orbitals. The La intraband 
transition is a difficult case to calculate with the unknown 
transition matrices involved. For simplicity and for the 
purpose of illustration, here we consider the part of the 
optical transition to the conduction bands with the H tc t 
character, and assume the transition matrix is non-zero 
only between the La-5<i and to its neighboring H tc t or- 
bitals. 



Rt 



b n (Rt)\n u s) 



(45) 
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Let |R t , s) be the H te t Is state at site R t , and b n (R t ) = 
(Rt, s\n) the amplitude of the state |Rj, s) in |n), we then 
have 



{vac \x\ n) — 



R La ,R t ,a 



l* a (Ria)&n(Rt)(RLa, d a \x\ R t , s), 



(46) 



where the sum over R t runs all the neighboring H te t of 
the six La sites. We calculate The a* (Rlq) and b n (R t ) 
by numerically solving a finite size cell (15 x 15 x 15oq) 
of LaH3 with a single H oct vacancy. The orbital overlap 
(d a |x|s) oc (d a \p x ) is a linear combination of the coupling 
parameters V p da and Vpdn and can be found in Harrison's 
book jil|. It is estimated the ratio Vpdv/Vpdn — —2.17. 
The optical conductivity (t(u>) can thus be calculated up 
to an overall constant. In the calculations, we use a rigid 
band approximation, where the conduction band spectra 
is unaffected by the vacancy. This approximation is jus- 
tified by comparing the result to the actual conduction 
band in the presence of the vacancy from the finite size 
calculation. Only minute deviation is observed in small 
size calculation. 

The calculated result (Lorentzian width 0.5 eV) is 
showed in the Fig. 16 which gives a fairly featureless 
curve. Except the small peak around 4 eV, the curve 
does not show any pronounced features. This is qualita- 
tively consistent with the experiments . 

There is a threshold ~ E c in a(w), which is the en- 
ergy cost for the vacancy state to transfer to the conduc- 
tion band. In our calculation, E c ~ 0.1 eV. This small 
value is due to the inaccurate tight-binding fitting to the 
LDA, resulting in lower energy at X-point in LaH 3 . If we 
subtract this, E c ~ 0.38 eV. In our estimate, a charged 
Coulomb attraction has not been added, so that the ac- 
tual threshold is larger. 

It is interesting to compare the optical properties of 
the hydrides vacancy state with those of the usual doped 
semiconductors studied by Thomas et al fljlf . First we 
note an important difference between the two systems. 
In a conventional semiconductor, the hydrogen-like Is 
groundstate has an electronic dipole active transition to 
the bound 2p state, leading to optical absorption peaks 
at about 30 meV in Si. In the metal hydrides, the ground 
state of a vacancy is a highly localized s-wave-like state 
(see Fig. |l^). 

However, the p- wave-like state is expected to be only 
weakly bound with much larger spatial distribution. This 
is because of the absence of short-range attraction as 
in the s-wave state. The long range Coulomb attrac- 
tion between the electron and the vacancy induces a 
bound p-wave state with a larger radius similar to the 
hydrogen-like ones. Therefore we expect only a weak 
dipole transition intensity from the very localized s-wave 
state (mainly induced by the short range bonding force) 
to the p-wave state. The absence of the short-range force 
in the p-wave vacancy state can be explained as follows. 



Consider a local p-wave-like state of La-5d orbits, sim- 
ilar to that in Fig. |l3| except for the symmetry that is 
changed from s- to p-wave-like. By symmetry, the p- 
wave-like state is decoupled to the Is electron state of 
the center hydrogen atom. Note that this property of 
the vacancy state is markedly different from that of the 
negative ion vacancy in an alkali halide crystal. The lat- 
ter is called a color center that has one excess electron 
bound at the vacancy Q . The potential created by the 
vacancy in alkali halides is a long-range Coulomb force, 
and the color center absorbs visible light in a dipole tran- 
sition to a bound excited state of the center. Our theory 
for the vacancies in hydrides predicts much weakness of 
such transitions. 



VII. CONCLUSION 

In this paper, we have examined the importance of the 
electron correlation in negatively charged ions (H _ ions), 
and identify the metal hydrides to be a strongly corre- 
lated system. We develop a many-body theory to de- 
scribe the H _ ion lattice in hydrides. We use lanthanum 
hydrides as a prototype for these compounds and show 
that LaH2 is a metal, and LaH3 is a band insulator where 
the electron correlation plays a crucial role. 

The electronic structure for metal hydrides is described 
by the large U -limit Anderson lattice model. The param- 
eters for the hydrides are determined by a combination 
of method including the microscopic calculations for the 
electron hopping integrals between hydrogen sites and 
parameter value obtained from LDA methods for La/H 
hybridization. The Anderson lattice model is then solved 
using a many-body technique — Gutzwiller method. The 
elementary entity when a H atom is removed from the in- 
sulating LaH3 is the octahedral hydrogen vacancy, which 
is highly localized electronically, due to the strong hy- 
bridization between the H-ls and its neighboring La-5cZ- 
e g orbitals. This explains why the metal-insulator tran- 
sition can occur at a large vacancy concentration ~ 25%. 

Our theory for the localized vacancies in hydrides pre- 
dicts a much weak optical transition to the hydrogen- 
like- 2p state in comparison with observed transitions in 
conventional semiconductors. 
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d(a.u.) 


4 


4.5 


6 


£i(a.u.) 


-0.25746 


-0.22687 


-0.16789 


a\ 


0.44324 


0.39654 


0.27863 


a 2 


0.08514 


0.05775 


0.01804 


0,3 


0.02785 


0.01376 


0.00146 


hi 


-0.14039 


-0.11154 


-0.05866 


b> 


-0.04145 


-0.02566 


-0.00748 


b 3 


-0.01609 


-0.00731 


-0.00123 


t 


-0.03611 


-0.02857 


-0.01346 



TABLE I. Values of integrals in Eq. (jlGj) and the hopping 
integrals t(d) for different inter-proton distances d. All dis- 
tances and energies are in atomic unit. 



Cubic size (ao) Crystal field reduction (eV) 

1 0.084674 

2 0.145003 

3 0.212051 

4 0.223285 

5 0.226712 

6 0.230330 



TABLE II. Change of ti, the electron hopping integral be- 
tween the n.n. tetrahedral H sites, due to crystal field effect of 
ions within a cube cell of various size. The change converges 
as the size of the cubic cell increases. a : lattice constant in 
LaH 3 . 
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Parameters 



atomic energy of La 5d ti g 
atomic energy of La 5d e g 
atomic energy of H to t 
atomic energy of H oct 
n.n. La-La 7r-bonding 
n.n. La-La cr-bonding 
n.n.n. La-La 7r-bonding 
n.n.n. La-La cr-bonding 
n.n. La-H oc t cr-bonding 
nearest La-Htet cr-bonding 
nearest H to t-Htct cr-bonding 
nearest H oct -Htct cr-bonding 



fitting values type of atom position (do) distance (cio) weighting 



C-d 



1.6 

e'd = 1.4 
et = -3.2 
e = -2.6 
V dd7r = 0.6 
Vu* = -1.2 

vU = -o-i 

V3,,,, = -0.2 

2 
1.3 

= -0.4 
= -0.79 



V'La — o i 



s d a 
La — 



TABLE III. Tight binding parameters extracted from the 
LDA results of Fig. ^ for LaH3 and LaH2. All the energies 
are in units of eV. The band structure of LaH3 from the tight 
binding model of these parameters is plotted in Fig. pi 



Parameters 


Renormalized values 


nearest La-H oct cr-bonding 
nearest La-H to t cr-bonding 
nearest H te t-H t ot cr-bonding 
nearest H oc t-H to t cr-bonding 


vj£-° = -1-004 

vfsr 1 = - 1 - 146 

V*~J = -0.228 
V£? = -0.309 



TABLE IV. Renormalized 
method for LaH3. 



parameters of Gutzwiller 



Htot 

Htet 

Htot 
Htot 

Htet 

Htot 

H oc t 
H ct 
Hoot 
H oc t 

La 
La 
La 
La 



)rr' 

\4> 4' 
(3 3 

(3 3 

V 4 ' 4 ' 
f5 I . 

V 4 ' 4 ' 
(5 3 
\A1 4i 

(- 1 

V 2 ' 2 ' 

(1,0, 

(M. 
(i,i, 
(|,o, 
(- 1 

V 2 ' 2 ' 



(|,o. 



-) 

0) 
0) 

I) 

0) 
0) 

I) 

0) 
0) 



0.4330 
0.8292 
1.0897 
1.2990 
1.2990 
1.4790 
0.7071 
1.0000 
1.2247 
1.4142 
0.5000 
0.8660 
1.1180 
1.5000 



0.0012 
0.0365 
0.0002 
0.0134 
0.0087 
0.0012 
0.1426 
0.0111 
0.0001 
0.0004 
0.5790 
0.0983 
0.0499 
0.0269 



TABLE V. Weighting factors of the vacancy state at differ- 
ent atomic sites. The third column shows the total weighting 
of atomic sites which share the same distance from the va- 
cancy site (origin). Only one position vector is shown in the 
second column for each distance while the others can be de- 
duced from symmetry. 



Vacancies separation (a ) 
y/2/2 (different s.l.) 
1 (same s.l.) 
3/2 (different s.l.) 
V2 (same s.l.) 
TO/2 (different s.l.) 



Hopping integral (eV) 



0.475 
-0.7 
0.125 
-0.175 
0.0625 



TABLE VI. Calculated hopping integrals of vacancy state 
for various vacancies separation. Inside the parentheses are 
the indication of the two vacancies belonging to the different 
or the same sublattices (s.l.) of the fee. 
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FIG. 1. Lattice structure (f.c.c) of LaEU. As x increases 
from 2 to 3 the H-atom content at octahedral sites increases 
from empty to full, and a shiny metal evolves to an insulator. 




FIG. 4. Ground state energies (E-, odd parity) and lowest 
energies of even parity state (E+) in a H^" ion as a function of 
the two proton distance d. The lower (upper) solid line shows 
the present result for the bonding energy E- (antibonding 
energy E+) of the HJ ion. Dashed line is the E- from Fis- 
cher- Hjalmars pH |. The lower dotted line is the E- estimated 
by Taylor and Harris Q|. The upper dotted line shows the 
ground state energy of a H2 molecule and a free electron for 
comparison. At d — * 00, the energy is the sum of a H~ ion 
and a H atom, shown with a square from the present calcu- 
lation, a triangle from Fisher- Hjalmars and the lower circle 
from Taylor and Harris. The upper circle is the energy of two 
independent H atoms and a free electron. 



FIG. 2. Illustration of Chandrasakhar's wavefunction 
Eq. (^) for H~, describing two electrons bound to a proton. 
The solid line between the two electrons represents the corre- 
lation term in Eq. (|l|). 
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CD 
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d (a.u.) 



6.0 



FIG. 3. Illustration of HJ ion, a system with two protons 
of distance d and three electrons. Shown at left is a H~ ion, 
described by the wavefunction of Eq. (|l|), and at right is a 
neutral H atom in the groundstate. 



FIG. 5. Electron hopping integral t(d) as a function of the 
distance d between two H ions. The hydrogen hopping in- 
tegral decreases rapidly with increasing distance between the 
two H-ions. The equilibrium separation between the two near- 
est Htct's, and between neighboring H tc t and H oct are 5.29 a.u. 
and 4.58 a.u., respectively. The solid line shown here is fitted 
to the three data points. 
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FIG. 7. Band structures of (a) LaH2 and (b) Lari3 from 
local density approximation calculation. The Fermi energy 
in LaH3 is between the top of the H bands (the lower 3 
bands) and the bottom of the La bands (upper 5 bands). 
T = (0,0,0), X= (2tt,0,0), W= (27r,7r,0), L= (7r,7r,7r), 
K= (3-7t/2, 3n/2, 0), in unit of l/a , are the high symmetry 
points in crystal momentum space. 
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FIG. 6. A cubic cell with the center of H^" as the origin 
is shown. Total charges of ions enclosed in the cell are neu- 
tral. Charges in unit of e on an edge, surface and corner are 
counted as a quarter, half and eighth of the original charges 
respectively. Only a few examples are given in the figure. 
Dotted lines are guides to the eyes. 
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FIG. 8. The fitted band structures from the tight binding 
model with parameters listed in Table for (a) LaH3 and 
(b) LaH 2 . 
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FIG. 9. The conduction and valence bands of LaH2 cal- 
culated from model Hamiltonian ( |l7|) by using Gutzwiller 
method. 
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FIG. 11. Band gap of LaH3 is observed to be a monotonic 
increasing function of lattice constant a (in unit of do). There- 
fore by exerting large enough pressure on the crystal, the 
band gap will be closed up and the originally insulating hy- 
dride will become conducting. This transition happens when 
a ~ 0.89a . 
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FIG. 12. 2-D illustration of an impurity state. Six neigh- 
boring La-5d-e 9 orbitals (only 4 of them are shown here) are 
pointing toward the extra H oct . 
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FIG. 10. Band structure of LaH3 of model Hamiltonian 
(jlTj) using Gutzwiller approximation. In (a) the crystal field 
is included in estimating the H-H hopping integrals. In (b) 
the crystal field effect is not included. The actual result is 
expected to be between (a) and (b). 
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+) 

FIG. 13. Diagrammatic illustration of the proposed H oc t 
vacancy state in LaH3. The center circle represents a H va- 
cancy, forming an n-type impurity center. The surround- 
ing orbits represent phases for the local s-like octahedral 
La-5d-e g . The vacancy state has spin- 1/2. 




FIG. 14. A 2-D illustration of vacancy state hopping. The 
dashed line with arrow represents the microscopic hopping 
process. Note that in 3-D, the nearest neighbor vacancy states 
have two common La atoms and have a shorter inter-vacancy 
distance. 
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FIG. 15. Density of states for (a) LaHh and (b) LaH3. Solid 
curves are our calculated results and dashed curves are the 
experimental data obtained from Peterman et al. Bq|. 




FIG. 16. Optical conductivity due to excitation of electrons 
from vacancy state to conduction band of LaH3. 
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